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^ ! Abstract 

T-H . Sound wave propagation in a relativistic perfect fluid with a non-homogeneous 

[ isentropic flow is studied in terms of acoustic geometry. The sound wave equation turns 

' out to be equivalent to the equation of motion for a massless scalar held propagating in 

■ a curved space-time geometry. The geometry is described by the acoustic metric tensor 

^ , that depends locally on the equation of state and the four-velocity of the fluid. For a 

Q\ ' relativistic supersonic flow in curved space-time the ergosphere and acoustic horizon 
may be defined in a way analogous to the non-relativistic case. A general-relativistic 

O , expression for the acoustic analog of surface gravity has been found. 

U '. 

1 Introduction 

> ■ 
• 1— I . 

^ ■ Since Unruh's discovery that a supersonic flow may cause Hawking radiation, the analogy 
. between the theory of supersonic flow and the black hole physics has been extensively dis- 
cussed 0, ^, ^, |, H, 0. All the works have considered non-relativistic flows in flat space-time. 
In a comprehensive review 0] Visser defined the notions of ergo-region, acoustic apparent 
horizon and acoustic event horizon, and investigated the properties of various acoustic ge- 
ometries. 

In this paper we study the acoustic geometry for relativistic fluids moving in curved 
space-time. Under extreme conditions of very high density and temperature, the velocity 
of the fluid may be comparable with the speed of light c. In such conditions, the speed of 
sound may also be close to c. In particular, the equation of state of an ideal ultrarelativistic 
gas yields the speed of sound as = c/v^S. Such conditions are physically realistic in 
astrophysics, early big-bang cosmology and relativistic collisions of elementary particles and 
ions. 

We restrict attention to perfect, isentropic fluids and, in addition, we impose that the flow 
be irrotational. In most of our considerations we make no assumption about the space-time 
metric. However, in the calculation of surface gravity we assume stationary space-time and 
steady flow. Therefore, in our discussion we make no distinction between acoustic apparent 
and event horizons. 
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We organize the paper as follows. In sections |^ and ^ we derive the relativistic acoustic 
wave equation in terms of an acoustic metric tensor. Section ^ describes the notions of 
acoustic horizon and ergo-region. In section |^ we give the calculation of surface gravity at 
the horizon and also study the acoustic Unruh effect and Hawking temperature. We draw 
our conclusions in section ^. In appendix ^ we outline the basic relativistic fluid dynamics 
and derive a few formulae needed for our consideration. 



2 Relativistic Acoustics 

We begin by deriving the sound wave equation for a perfect, irrotational gravitating fluid. 
We denote by , p, p, n and s the velocity, pressure, energy density, particle number 
density and entropy density of the fluid. The basic equations of relativistic fluid dynamics 
(for details see appendix 0) needed for our discussion are the continuity equation 

d^iy^nu'') = (1) 

and the equation of potential flow 

wUf, = -d^(t) , (2) 
where the specific enthalpy w is defined as 

^ = -^. (3) 

Following the standard procedure [0, |^ we linearize equations (|l]) and by introducing 

w w + 6w, n ^ n + 6n, 

^ + ^ + 50, (4) 

where the quantities 6w, 6n, 5u^ and 50 are small acoustic disturbances around some average 
bulk motion represented by w, n, and 0. For notational simplicity, in our equations we 
shall use ip instead of 50. Fluctuations of the metric due to acoustic disturbances will be 
neglected. Note that the normalization ([T^ ) implies 

g^uu^du" = 0. (5) 

Substituting into (Q) and retaining only the terms of the same order of smallness, we find 

5w = —u'^d^ip (6) 

and 

wSu" = -g^'^d^ifi + u^'u^d^ifi. (7) 
Similarly, linearization of the continuity equation (|l|) gives 

d^,{^/^6nu^) + d^iy^ndu^") = 0. (8) 

Replacing 6n by 

Sn={^r'6w, (9) 
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and substituting (|) and (|^) into 



n 
w 



-9 



we obtain the wave equation 
^ n dw 



1 



'wdn 



d^Lp = 0. 



(10) 



This equation describes the propagation of the hnearized field ip. The space-time dependence 
of the background quantities n, w and = — w~^c?^0 is constrained by the equations of 
motion for a perfect, isentropic and irrotational gravitating fiuid. 

In the simplest case of a homogeneous fiow in fiat space-time, equation (|1^) becomes the 
free wave equation 

(d^ - c^A)^ = 0, (11) 



where the quantity Cg is the sound wave velocity given by 



c = — 



n dw 
w dn 



(12) 



s/n 



The subscript s/n denotes that the derivative is taken while keeping the specific entropy s/n 
constant, i.e. for an isentropic process. By making use of equation (|80D from appendix ^ 
it may be easily verified that this definition is equivalent to the standard definition of the 
speed of sound [H 

' dp 

Introducing the symmetric tensor 



(13) 



s/n 



r = 

w 



(1 



we rewrite the wave equation (|T0|) as 



(14) 



(15) 



This form of the acoustic wave equation will be used to construct the acoustic metric. 

At this stage it is instructive to calculate the tensor f^iy in the fiat-metric non-relativistic 
limit. To do this, we reinstate the speed of light c and make the replacements 
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w - 



mc 



7]^ 

vjc 



n 



.(l,t;Vc), 



(16) 



in (P^). Retaining only the leading terms in 1/c, we find 

Pnr 



/ 



^NR 2 2 2 



(17) 



This 4x4 matrix, apart from the overall constant factor l/(mc)^ and the minus sign due to 
the signature convention, precisely equals the corresponding non-relativistic tensor discussed 
by Unruh H] and Visser P. 
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3 Acoustic metric 



We now show that equation (|T5D with (|1J) can be put in the form of the d'Alembertian 
equation of motion for a massless scalar field propagating in a (3+l)-dimensional Lorentzian 
geometry described by the acoustic metric tensor G^^. We proceed in a way analogous to 

i- 

We first introduce the acoustic metric tensor G^j^p and its inverse G^'^ such that 



where 

These equations imply 



G = detif^"). 
det(G^^) = det(G"^'^))-^ = G. 



(19) 
(20) 



It follows from equations (|1^) and (|T^ that the metric G^i, must be of the form 



G,, 



'i^u - k[g^p - (1 - cJm^Mj.], (21) 
where is a conformal factor which will be fixed later. Indeed, one may easily check that 



using this expression and (plSf ) one obtains 



G^'^G^u — k — — 6^ . 

Wy/—G 

To calculate the determinant of G^u, it is convenient to use comoving coordinates 
the general comoving coordinate system, the four-velocity vector is given by 



U,, 



and hence 



G = k^det L,-(l-a^) 



(22) 
1- In 

(23) 
(24) 



This expression may be simplified using standard algebraic properties of a determinant. We 
find 

G = k^clg. (25) 
This equation combined with (p2|) fixes the factor k as 

n 



k 



wc. 



Thus, we have shown that the acoustic wave equation may be written in the form 

^ ^,i^/^Gn^uV = 0, 



with the acoustic metric tensor given by (^Tj) and its inverse by 



(1- 



(26) 



(27) 



(28) 
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4 Ergosphere and acoustic horizon 



In analogy to the Kerr black hole we may now define the ergo-region as a region in space-time 
where the stationary Killing vector becomes space-like. The magnitude ||^|| of may be 
calculated in terms of the three velocity v defined in appendix 0. Using equation (|2T|) (and 



(M)- (p6|) from appendix ^ we find 

iieip = G,.ec = keu^ - (1 - c^)f ] = ke^,i\c'^ - v^). (29) 

This becomes negative when the magnitude of the fiow velocity v exceeds the speed of sound 
Cg. Thus, as in non-relativistic acoustics |0, any region of supersonic fiow is an ergo region. 
The boundary of the ergo-region is a hypersurface defined by the equation 

v"^ -cl = 0. (30) 

Since the stationary Killing field becomes null on , this hypersurface is called stationary 
limit surface fl^. The hypersurface may, in general, be quite complicated, in particular 



if the fiuid equation of state, and hence the speed of sound, is space-time dependent. For 
simplicity, in the following discussion we assume Cg = const throughout the fiuid. 



To define acoustic horizon we use the concept of wave velocity of a hypersurface [|Tl|]. Let 
{Sq} denote a set of hypersurfaces defined by 

-a^ = 0, (31) 

where a is a constant, < a < c. The stationary limit surface is contained in {S^}. 



Equation (|3lD together with the equation 

xq = const (32) 

determines a two-dimensional surface aa- The normal to each is given by the vector field 

oc df.v'^. (33) 

Now we define 

v± = ~ , = (34) 

as the wave velocity of two-dimensional surfaces (Jq as measured by the observer u, whose 
world-line is such that the four-velocity evaluated at a point of is tangent to the world- 
line of u. Clearly, v'j_ < if and only if is time-like or null, i.e. if and only if n^n^ < 0. 
If Sa is time-like, may be normalized as n^n^ = — 1. In this case, we decompose as 

u^" = v±j±n'' + 7xL'', (35) 

where 

7± = -A=^^ (36) 
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and 

L'' = —{g>^- + n^nnu, (37) 

7± 

is a time-like unit vector which represents a displacement in along the projection of u^^. We 
can regard (|35|) as a decomposition of the fluid flow into normal and tangential components 
with respect to Sq. The tangential three- velocity f ^ may be found by decomposing the vector 
in a way similar to (|85D from appendix 0: 

L'^ = Tii^'' + ig'"" - tn')L,, (38) 

where 



711 = t,L^ = (39) 

and 



vi = —V ; i = 1, 2, 3 . (40) 
For steady flow, t^n^ = 0, and hence equation (P^D, together with (RBI) and (B^), yields 



7 = 7||7±- (41) 
Now we define the acoustic horizon as a hypersurface 7i defined by the equation 

^ c2 = 0, (42) 



i.e. a hypersurface the wave velocity of which equals the speed of sound at every point. The 
acoustic horizon Ti and the stationary limit surface in general do not coincide. Equation 



(PI) states that 7i is inside the ergo-region and overlaps with at those points where 
= t^, i.e. at the points where the three velocity is perpendicular to the-twosurface ac^. 
Let us illustrate our general considerations by an example. Consider a two-dimensional 
axisymmetric flow in a stationary axisymmetric space-time. Physically, this may be a model 
for the equatorial slice of a rotating star, or for a draining flow in an axially symmetric 
bathtub. In the coordinates t, (p, r, the general form of a stationary axisymmetric metric is 
given by 

rfs^ = -V{r){dt - W{r)d(j)y - V{r)-^r^d(j)'^ - nir^dr^, (43) 

where V, W and Q are arbitrary functions of r approaching V,Q ^ 1 and W ^ as r 
tends to infinity. We further assume that these functions are positive and non-singular in 
the region occupied by the fluid. According to equation (0), the fluid velocity takes the 
form 

u^ = ^{V-'^^-Wv^;v^,Vr) , (44) 

where 

= V-^r%l + (45) 

with f0 and Vr being the r-dependent azimuthal and radial components, respectively. The 
surfaces (Tq defined by equation ( PT| ) with ( [5^ ) are concentric circles with the centre at the 
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origin. The unit vectors n^^ and appearing in the decomposition of in (|35| ) may now 
be easily constructed using (0) and (|37D, respectively. We find 



= (0; 0, 1) , (46) 
L'' = 7\\{V/'-Wv^;v^,0) , (47) 

where 

vj = IZ-V^J. (48) 
From equation (|3lD, the wave velocity f j_ of hypersurfaces cTq is 



f_L = 'y\\flVr ■ (49) 

The equation of continuity and the condition that the fiow be irrotational constrain the 
radial dependence of v± and v\\. From equation (0) we find 

n'y±v±r = ci (50) 

and similarly, as a consequence of equation (p2D, we obtain 



m;7_l7h ( ^^^11 — 1 = C2, (51) 

where Ci and C2 are arbitrary constants. Of course, these expressions are valid as long as 
v'j_ and vjj do not exceed the speed of light squared = 1. If we further assume that the 
particle number n and the speed of sound Cg are constant throughout the fiuid, equation ( ^0|) 
yields an explicit expression for the horizon radius 



1-d 



Ci 



r-H = ^ - — . (52) 

C, n 



The radius of the stationary limit surface can be found by solving the equation 

1/2 

(53) 



1 -c 
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Evidently, the acoustic horizon coincides with the stationary limit surface, i.e. the boundary 
of the ergosphere, if both the azimuthal velocity and the off-diagonal metric term W{r) 
vanish. It is remarkable that the horizon radius th does not depend on metric, whereas the 
radius of the ergosphere does. 



5 Acoustic surface gravity and Hawking temperature 

Having established the definition of the acoustic horizon, we proceed to calculate surface 
gravity at the horizon, restricting attention to stationary acoustic geometries. If the fiow is 
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steady, and the metric g^i, is stationary, there exists a system of coordinates in which the 
components of G^i, are time independent and the acoustic geometry is said to be stationary. 
We shall define the surface gravity in terms of the Killing field that is null on the 



horizon. The surface gravity n may then be calculated using [|12 

G^'^d.WxW' = -2^X' (54) 

evaluated on the horizon, where 

llxir = G'„;3xV. (55) 

Since = on the horizon, the derivative c^^HxlP is normal to the horizon, so that on 

the horizon we have 

G'^^r^,|-||x|P = 2«:x^ (56) 

where d/dn = n^d^ denotes the normal derivative on the horizon. The definition of surface 
gravity by equation is conformally invariant [0. Therefore, in the calculations that 
follow we drop the conformal factor k in G^jj. 

First, consider the simplest case when the stationary limit surface and the acoustic 
horizon Ti coincide. Then, x^ is equal to the stationary Killing field which, according to 
equation (|29|) , is null on the horizon. A straightforward calculation gives 

1 1 

G^^n^ = n^ u^" = ^= (57) 

on the horizon. Equation ( ^6]) applied to together with ( ^7|) and (^91) yields 



where it is understood that the derivative and the quantity ^''^i^ are to be taken at the 
horizon. The corresponding Hawking temperature Th = k/{2'k) represents the temperature 
as measured by an observer near infinity. The locally measured temperature follows the 
Tolman law, i.e. 



so that the local temperature approaches 

1 d 



c? dn 



(60) 



as the acoustic horizon Ti is approached. Thus, in this limit, equation (|59[) with (^) corre- 
sponds to the flat space-time Unruh effect. This behaviour is analogous to the Schwarzschild 



black hole [0] except that our local temperature does not diverge as one approaches Ti. 

Next, we consider the situation when Ti and do not coincide. In this case, the 
calculation of surface gravity may, in general, be quite non-trivial. As in non-relativistic 
acoustics 10, there is no particular reason to expect the acoustic horizon to be in general a 
Killing horizon. However, the matter may greatly simplify if we assume a certain symmetry, 
such that an analogy with the familiar stationary axisymmetric black hole may be drawn. 

We consider the acoustic geometry that satisfies the following assumptions: 
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1. The acoustic metric is stationary. 



2. The flow is symmetric under the transformation that generates a displacement on the 
horizon along the projection of the flow velocity. 

3. The metric Qf^i, is invariant to the above transformation. 

The last two statements are equivalent to saying that a displacement along the projection 
of the flow velocity on the horizon is an isometry. The axisymmetric geometry considered in 
section ^ is an example that fulfils the above requirements. 

Suppose we have identified the acoustic horizon as a hypersurface TC defined by equation 
(^). Consider the vector field defined by (^). Its magnitude is given by 

\\L\\' ^ G,,L^L^ = [1 - (1 - clWA = llicl - vl) (61) 

Hence, the vector L'^ is null on the horizon. We now construct a Killing vector x'* in the 
form 

x'' = e+ujr, (62) 

where the Killing vector ip^ is the generator of the isometry group of displacements on the 
horizon, and the constant u is chosen such that becomes parallel to L'^ and hence null on 
the horizon. The quantity u is analogous to the horizon angular velocity of the axisymmetric 
black hole [0 . 

We use a special coordinate system in which the Killing vector ip has the components 

r = (63) 

where the subscript ip denotes the spatial coordinate x'^ associated to the above-mentioned 
isometry. The metric components g^^, and the flow velocity field are, by assumptions 2 and 
3, independent of x^. In this coordinate system, the vector takes the form 

L^ = ^J^-'-^;v,A0] . (64) 



'9oo 9oo 

Requiring this vector to be parallel to the horizon, we set 



= VrT.L^, (65) 

with 

vrt;;=7,r'f-L-*^)". (66) 

" \V9oo 900 J 



This equation combined with (|H2|) determines the "angular velocity" of the horizon 



.=i^-'-^y\,. (67) 

fi'oo 9oo / 



The calculation of the surface gravity is now straightforward. The magnitude of in 
the vicinity of the horizon is given from equation (^) by 

\\x\\' = ll{cl-vl)x''Xu. (68) 
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Using the definition of (p7| ) and equation (0), we find 



1 I 

Gf^^n^ = 71'' = (69) 



on tlie liorizon. Equation (p6D, togetlier with and (pOf), yields 

where the values of the derivative and the quantity y^Xv are to be taken at the horizon. 
The non-relativistic limit of this equation in flat space-time coincides with the expression for 
surface gravity derived by Visser 0]. 



6 Concluding Remarks 

We have studied the relativistic acoustics in curved space-time in terms of an acoustic metric 
tensor. The acoustic metric tensor involves a non-trivial explicit space-time dependence. 
Apart from that, all the features of the non-relativistic acoustic geometry in flat space-time 
0, 1^ have their relativistic counterparts. We have shown that the propagation of sound 
in a general-relativistic perfect fluid with an irrotational flow may be described by a scalar 
d'Alembert equation in a curved acoustic geometry. The acoustic metric tensor is a matrix 
that depends on the fluid equation of state and involves two tensors: the background space- 
time metric tensor g^^^y and the product of the fluid four- velocities u^Uy. We have then shown 
that if in a non-homogeneous flow there exists a supersonic region called ergo-region, then 
an acoustic event horizon forms. We have discussed an example of axisymmetric geometry 
where the acoustic horizon does not coincide with the boundary of the ergo-region. We have 
calculated the surface gravity k, using a conformally invariant definition which involves 
a properly defined Killing vector that is null on the horizon. An acoustic horizon emits 
Hawking radiation of thermal phonons [Q] at the temperature Th = k/{27t) as measured by 
an observer near infinity. 

We believe that the effects we have discussed may be of physical interest in all those 
phenomena that involve relativistic fluids under extreme conditions. This may be the case 
in astrophysics, early cosmology and ultrarelativistic heavy-ion collisions. 

A Basic Fluid Dynamics 

Consider a perfect gravitating relativistic fluid. We denote by , p, p, n and s the velocity, 
pressure, energy density, particle number density and entropy density of the fluid. The 
energy-momentum tensor of a perfect fluid is given by 

Tf,u = {p + p)Uf,Uu - pg^y, (71) 

where g^^ is the metric tensor with the Lorentzian signature (H ). Hence, in this 

convention, we have 

u^u^ = g^yuf^u" = 1. (72) 
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Our starting point is the continuity equation 

{nu^),, = -^d^{^ nu^') = 0, (73) 
and the energy-momentum conservation 

T^r^^ = 0. (74) 
This equation apphed to (|7l| ) yields the relativistic Euler's equation |^ 

{p + p)u''u^-u + d^p + u^u^duP = 0. (75) 

Euler's equation may be further simplified if one restricts consideration to an isentropic fiow. 
A fiow is said to be isentropic when the specific entropy s/n is constant, i.e. when 

d,{^) = 0. (76) 

A fiow may in general have a non-vanishing vorticity uj^u defined as 

uJ^^^ = hP^hyip.^], (77) 

where 

K = 6i:- u^u, (78) 

is the projection operator which projects an arbitrary vector in space-time into its component 
in the subspace orthogonal to u^. A fiow with vanishing vorticity, i. e. when 

= 0, (79) 

is said to be irrotational. In the following we assume that the fiow is isentropic and irrota- 
tional. 

As a consequence of equation (|76| ) and the thermodynamic identity 

s 1 

dw = Td(-) + -dp, (80) 
n n 

with w = {p + p)/n being the specific enthalpy, equation (|75D simplifies to 

u''{wup),,-d^w = 0. (81) 
Furthermore, for an isentropic irrotational fiow, equation ( |79D implies ||TT 



{wUf,).y - {wu^).f, = 0. (82) 
In this case, we may introduce a scalar function such that 

wu^ = -(9^0, (83) 
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where the minus sign is chosen for convenience. It may be easily seen that the quantity wu^ 
in the form (|83D satisfies equation (|8l|) . Solutions of this form are the relativistic analogue 
of potential flow in non- relativistic fluid dynamics [Q. 

It is convenient to parameterize the components of the fluid four- velocity in terms of three- 
velocity components. To do this, we use the projection operator g^i, — t^t^, which projects a 
vector into the subspace orthogonal to the time translation Killing vector = (1; 0) , where 



is the unit vector 



We split up the vector in two parts: one parallel and the other orthogonal to t^: 
where 

7 = f'u^ ■ 

From (|85|) with (^) we find 

uo = l^/goo- 

Equations (pS]), and the normalization of imply 

-fiju'u^ =7^-1, 
where 'jij is the induced three-dimensional spatial metric: 

lij = --Qio'i «,J = 1,2,3. 

5'oo 

Now, it is natural to introduce the three- velocity w*, so that 

= 'yv\ 

with its covariant components Vi and the magnitude squareed f ^ given by 
It follows from (|88|), and that 

1 
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92) 



Since u'^ and are time-like unit vectors, a consequence of (^) is that 7 > 1 and hence 
0<v^<l. 

Equations (Q) - (pl| ) may now be used to calculate other covariant and contravariant 
components of the fluid four- velocity. We find 
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